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OVER  AEROFOILS  WITH  NOVEL  REAR  SECTIONS 

by 

P.  Ashill 


SUMMARY 

C- 

A  combined  theoretical  and  experimental  investigation  into  transonic  flows 
over  aerofoils  of  advanced  design  is  described.  The  experiments  have  been'^ 
performed  at  high  subsonic  speeds  and  over  a  wide  range  of  Reynolds  number  up  to 
20  million  on  a  number  of  aerofoils  with  rear  pressure  distributions  of  differing 
form  and  severity.  Three  families  of  aerofoils .have  been  studied,  all  of  which 
are  of  14%  thickness  and  have  a  high  degree  of  rear  camber.  Data  for  aerofoil 
pressure  distribut ions  and  overall  forces  are  compared  with  predictions  by  a 
calculation  method  based  on  the  viscous-inviscid  interaction  concept  and 
including  allowance  in  the  modelling  of  the  turbulent  shear  layers  for  effects 
which  become  important  as  separation  is  approached.  Predictions  of  pressure 
distributions  by  this  method  are  shown  to  be  in  generally  good  agreement  with 
experiment  for  the  various  sections  tested,  and  it  is  concluded  that  the  method 
provides  a  reliable  basis  for  designing  aerofoil  sections  with  advanced  features. 
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I  Introduction 

CompijfcaC  lonal  Fluid  Dynamlna  (CFD) 
methods  Tor  transonlo  riows  over  aerofoils 
have  been  shown  to  yield  physiaally 
r-ali.iti.;  solutions  in  a  number  of  special 
onaes  generally  with  fully  attached  flow. 
However,  these  methods  treat  the  shear 
layers  approximately  and  hence  further 
assessment  Is  needed  before  the  concepts 
Involved  can  he  extended  with  confidence 
to  more  complex  flows.  The  need  for  such 
an  asses.sment  Is  particularly  acute  for 
flows  at  high  subsonic  speed  over  aero¬ 
foils  of  modern  design  having  relatively 
large  rear  camber.  For  such  flows,  the 
adverse  pressure  gradients  over  the  rear 
of  the  section  can  be  particularly  severe, 
leading  to  a  rapid  growth  In  boundary- 
layer  thickness  as  Che  trailing  edge  is 
approached  and  possibly  also  to  separ¬ 
ation.  Thus  flows  of  this  type,  which  are 
of  great  practical  Importance,  pose  a  con¬ 
siderable  challenge  to  the  accuracy  of  the 
predictions  by  CFD  methods. 

This  paper  describes  a  combined  theor¬ 
etical  and  experimental  investigation 
aimed  at  Improving  the  modelling  In  a  CFD 
metiiod  of  flows  approaching  separation. 

I'he  experiments  were  performed  at  high 
subsonic  speeds  In  the  8ft  x  8ft 
Pressurised  Wind  Tunnel  at  RAE  Bedford  on 
three  families  of  aerofoils  of  ItX 
thickness  and  of  modern  design,  some  of 
which  are  novel  In  character. 


Until  fairly  recently,  the  method  most 
favoured  in  fJK  for  calculating  transonic 
flows  over  aerofoils  was  the  viscous  ver¬ 
sion  of  the  darabedlan  and  Korn  program 
(VGK)i.2.  Based  on  the  vlscous-lnvlsold 
Interaction  concept,  VGK  j.ses  a  full- 
potential  algorithm  for  the  Invlscld  flow 
and  Integral  methods  to  calculate  the 
shear  layers.  VGK  gives  reasonable  pre¬ 
dictions  of  aerofoil  pressure  distri¬ 
butions  and  drag  for  attached  flows  but, 
where  flow  separation  la  approached,  the 
method  Is  less  reliable.  The  Inadequacy 
of  the  method  for  flows  of  this  tvne  wems 
In  part  from  the  use  of  the  concepts  of 
first-order  boundary-layer  theory. 

More  recent  developments  have  Included 
methods  for  solving  the  Euler  equations’.^ 
In  the  Invlscld  part  of  the  flow  and 
others  yielding  solutions  of  the  Reynolds- 
averaged  Mavler-Stokes  equations’.’.  The 
former  class  of  methods  retains  the 
viscous-lnvlsc Id  Interaction  procedure  and 
thus  needs  allowance  for  effects  which 
become  Important  In  regions  approaching 
separation;  the  latter  class  Is  ‘  ss 
restricted  In  this  sense  but  currently 
relies  on  primitive  turbulence  models’. 
Thus  methods  of  this  type  are  not,  at 
present,  suitable  for  routine  application. 

In  order  to  remedy  the  defects  In  ’.’GK 
mentioned  above.  Improvements  have  been 
made  to  the  method  to  allow  for  ’higher- 
order’  effects  and  other  Influences  which 
become  important  In  separated-f low 
regions.  A  first  attempt  to  Include  .these 
effects  (to  give  a  code  known  as  AVGK)  was 
not  entirely  successful  because  of  the 
slowness  of  the  numerical  convergence  of 
the  'direct'  method  (used  to  couple  the 
Invlscld  and  viscous  flow  solutions)  as 
separation  Is  approached.  East  et  al’ 
showed  that.  In  such  circumstances,  an 
'inverse'  scheme  for  calculating  the  shear 
layers  Is  preferable.  In  this  scheme, 
boundary-layer  displacement  thickness  (or 
’wall'  transpiration  velocity)  Is  assumed 
known  and  the  wall  streamwls.e-veloclty 
distribution  of  the  Invlscld  flow  Is 
calculated;  this  distribution  Is  then  used 
In  the  invlscld  flow  either  directly  as 
the  wall  boundary-condition  or  indirectly 
via  an  intermediate  expression  to  define 
the  wall  transpiration  velocity.  The 
latter  'semi-inverse'  procedure  Is 
employed  In  an  Improved  version  of  AVGK 
known  as  BVGK'* ;  thus  AVGK  Is  superseded  by 
BVGK  and  so  no  further  mention  is  made  of 
AVGK  in  this  paper. 

The  remainder  of  this  paper  starts  with 
a  description  in  section  2  of  the  aero¬ 
foils  studied,  follows  this  with  a 
discussion  of  the  experiment,  including 


*:ne  ?:esc  techniques,  In  section  3,  and 
ionclnues  hy  iescrlblng  the  main  features 
cf  BVGK  In  section  Finally,  predic¬ 

tions  by  7GK  and  3VGK  of  aerofoil  pressure 
1 1s  c rlhiic  Ions  and  drag  are  compared  with 
those  measurement  In  section  5* 

2  Aerofoil  Design 

'\3  noted  In  the  Introduction,  three 
families  of  aerofoils  have  been  tested, 
and  In  this  section  characteristics  of 
eacn  family  are  briefly  summarised.  Each 
family  Is  11 s t IngM ished  by  the  pressure 
ilstrlbucion  or  shape  of  the  rear  half  of 
the  upper  surface.  In  other  respects  the 
aerofoils,  whicn  have  a  common  section 
'^head  of  65X  chord,  are  similar.  Details 
of  section  shape  and  rear  pressure  dlstri- 
oution  calculated  by  BVGK  for  selected 
aerofoils  are  given  In  Pig  1. 


«  Qtar  prtssuft 
disiribuMon 


■  n  1  0  735.  «  .  !0  .  10*.  Cl  «0  1 


lift  coefficient  -  0.6,  a  free-atreain 
Mach  number  M,  =  0,735,  a  chore  Reynolds 
number  n  »  2C  «  10°  and  with  transition  at 
5t  chord  on  both  surfaces,  separation  Is 
Just  avoided  near  the  trailing  “dge  on  the 
upper  surface.  RAE  5230  is  the  most 
extreme  of  the  three,  having  been  leslgned 
to  have  boundary-layer  characteristics 
near  the  trailing  edge  on  the  upper  sur¬ 
face  at  R  w  20  «  10°  similar  to  those  of 
RAE  5225  at  R  =  6  -  10°. 

2.2  Relaxing 

The  name  relaxing  Is  used  for  these  sec¬ 
tions  because  the  pressure  gradient  cn  the 
upper  surface  decreases  or  'relaxes'  with 
streamwlse  distance  between  751  and  90t 
chord.  This  type  of  section  has  a 
possible  advantage  over  more  conventional 
sections  for  wings  designed  for  high-speed 
cruise,  since  the  upper-surface  shape 
between  55X  and  75*  chord  is  such  that 
shocks  forming  in  this  region  are  rela¬ 
tively  weak  at  a  given  angle  of  incidence. 
Of  the  two  relaxing  sections,  RAE  5236  Is 
the  more  severe  In  terms  of  proximity  to 
separation  but  is  less  prone  to  separation 
than  all  the  convex  sections.  Despite 
this,  HAS  5235  has  a  higher  drag  at  a 
given  lift  at  *  0.735  than  does  the 
least  extreme  convex  section  RAE  5225,  as 
is  shown  later. 

2.3  Two  Part 

The  pressure  distribution  on  the  upper 
surface  of  the  solitary  two-part  section, 
RAE  523A,  Is  characterised  by  a  region  of 
modest,  adverse  pressure-gradient  between 
60*  and  90*  chord  followed  downstream  by  a 
pronounced  increase  in  pressure  gradient. 

A  novel  feature  of  this  section,  which  has 
a  }*  thick  base.  Is  that  it  has  a  rela¬ 
tively  large  cross-sectional  area  aft  of 
55*  chord,  offering  possible  structural 
advantages  and  some  freedom  to  design  a 
flap  with  a  generous  nose  radius. 

3  Wind  Tunnel  Tests 


1  Familic.  o7  aerofoils  sfudied 
2.1  Convex 

The  convex  family  of  sections  is  so 
named  because  the  pressure  gradient  on  the 
upper  surface  increases  monotonloally  with 
streamwlse  distance  between  about  60*  and 
80*  chord.  Six  sections  have  been 
Investigated  In  this  family,  three  with 
sharp  trailing  edges,  RAE  5225,  RAE  5229 
and  RAE  5230,  and  the  remainder  with  blunt 
bases.  However,  only  those  sections  with 
sharp  trailing  edges  are  considered  In 
thl.s  paper. 

or  the  three  sharp  tralllng-edge  aero¬ 
foils,  RAE  5225  has  the  least  tendency 
towards  separation  on  the  upper  surface; 
this  section  Is  designed  so  that,  at  a 


3.1  Model  and  Measurement  Techniques 

The  model  was  of  0.635  m  chord  for  all 
aerofoil  sections  and  spanned  the  tunnel 
working  section,  giving  an  aspect  ratio  of 
3.89  and  a  chord  to  height  ratio  of  0.26. 
Special  efforts  were  made  to  seal  the  two 
spanwlse  extremities. 

Boundary-layer  transition  was  fixed  by 
means  of  the  alr-lnjeotlon  technique  in 
which  air  Is  bled  into  the  boundary  layer 
thr-ugh  a  row  of  small  holes,  drilled  nor¬ 
mal  to  the  surface,  at  a  rate  Just  suf¬ 
ficient  to  fix  transition.  Details  of  the 
system  are  given  in  Pig  2  which  also  shows 
the  construction  of  the  model;  further 
Information  is  provided  In  Ref  10  where  It 
is  shown  that,  above  a  critical  mass  flow, 
transition  Is  effectively  fixed  at  the  row 
of  air  holes  at  5*  chord  on  both  surfaces, 
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ind  “vl’ieno"  Is  presented  suKgestlni;  *-,hat 
voe  npur'iouG  drag  sssocla^'ed  with  air 
Injection  Is  negligible. 

Surface  static  pressures  were  Tieasured 
at  orifices  of  0.5  mm  diameter  and  drilled 
normal  to  the  surface  both  on  the  model 
centre  line  and  at  a  limited  number  of 
positions  on  lines  about  one  chord  either 
side  of  the  centre  line.  The  off-centre 
holes  were  provided  as  a  check  on  the 
spanwlse  uniformity  of  the  flow  but  In 
this  paper  only  pressures  measured  at  the 
centre-line  are  considered. 

Section  lift  and  Ditching  moment  were 
determined  by  appropriate  Integrations  of 
the  static  pressures  around  the  section 
contour.  Section  drag  was  determined  from 
th“  momentum  deficit  In  the  wake  far 
downstream;  this  deficit  was  Inferred  from 
measurements  of  total  and  static  pressures 
made  with  a  rake  of  pitot  and  static  tubes 
located  at  the  vertical  plane  of  symmetry 
about  two  chords  downstream  of  the 
trailing  edge. 

Details  of  the  corrections  applied  to 
the  data  for  static-hole  error,  tunnel- 
wall  constraint  and  model  static  aero- 
elasclc  distortion  are  given  In  Ref  11. 

All  the  data  presented  In  this  paper  are 
corrected  for  these  effects.  In  particu¬ 
lar,  a  simple  correction  for  blockage  Is 
applied  to  free-stream  Mach  number  and 
static  pressure  on  the  basis  that  the 
blockage  Increment  tn  Mach  number  varies 
slowly  along  the  aerofoil  chord.  On  the 
other  hand,  a  wall  constraint  correction 
to  model  Incidence  alone  Is  not  Justified 
because  of  significant  variations  In  wall- 
induced  upwash  In  the  region  of  the  model. 
Therefore,  In  the  comparisons  between 
calculation  and  measurement  to  be 
discussed  In  section  5,  an  allowance  is 
made  for  this  variation  by  a  correction  to 
the  aerofoil  camher  In  the  calculation, 
using  the  analogy  between  flow  curvature 
and  aerofoil  camber  suggested  by 
linearised  aerofoil  theory '2, 

The  correction  for  model  aeroelastlc 
distortion  arises  because  the  model  was 
free  to  rotate  In  a  bearing  at  one  span- 
wise  extremity  while  being  constrained  at 
the  other  end  by  the  Incidence  adjustment 
device.  This  arrangement  avoided  the 
complication  of  a  following  mechanism  but 
meant  that  the  model  twisted  under 


aerodyharalc  load.  The  correction  to 
centre-line  Incidence  was  deduced  from  a 
prior  static-calibration  and  from  values 
of  aerodynamic  pltchlng-moment  Inferred 
from  the  static-pressure  measurements. 

The  same  calibration  Indicated  that  the 
twist  over  the  central  50t  span  was  about 
half  the  change  In  centre-line  incidence. 

Limited  checks  were  made  of  the  spanwlse 
variation  of  section  drag  coefficient  by 
rotating  the  wake  rake  about  Its  axis 
thr"ug.h  the  wake  for  a  number  of  shock- 
free  flows.  Differences  between  the 
centre-line  value  and  values  at  between 
0.3  and  0.5  chords  either  side  of  the 
centre-line  were  found  to  be  small,  typi¬ 
cally  from  0.0001  to  D.0002. 

An  assessment  of  the  effect  of  various 
sources  of  error  suggests  that  static 
pressure  and  drag  coefficients  are 
accurate  to  within  i0.002  and  lO.OOOl, 
respectively.  A  high  stanuard  of  repeat¬ 
ability  was  achieved  In  the  measurement  of 
pressure  both  between  and  within  test 
series,  static  pressure  and  drag  coef¬ 
ficients  being  repeatable  to  within  tO.OOl 
and  tO.OOOl. 

3.2  Test  Conditions 

The  measurements  were  made  for  free- 
streara  Mach  numbers  within  the  range  O.p 
to  0.75  and  for  chord  Reynolds  numbers 
from  6  «  10°  to  20  «  10°  for  all  sections 
except  RAE  5230  for  which  the  maximum 
Reynolds  number  was  17.7  «  10°.  All  the 
tests  described  here  were  made  with  tran¬ 
sition  fixed,  established  by  varying  alr- 
Injectlon  mass  flow  until  there  was  no 
further  change  In  drag  with  mass  flow. 

A  Calculation  Method 

BVGK  contains  the  same  basic  procedures 
as  VGK,  Including  the  lig-entralnment 
method  1 5  for  calculating  the  turbulent 
shear-layers,  but  embodies  changes  to  the 
modelling  of  the  shear  layers  which  become 
increasingly  Important  as  separation  Is 
approached  and  uses  a  revised  method  to 
couple  the  viscous  and  Invlscld  flow 
solutions . 

It .  1  Improvements  to  the  yiodelllng  of 
Viscous  Effects 

The  changes  to  the  modelling  of  the 
shear  layers  In  BVGK  are: 

(a)  Allowance  for  'higher-order'  effects 
In  the  equations  foi  Integrated  shear 
layer  mass  and  momentum  deficit  due 
to  flow  curvature  and,  where 
appropriate,  second-order  Reynolds 
stresses  . 

(b1  A  revised  shape-parameter  relation¬ 
ship  which  Is  more  suitable  for  flows 
with  separation. 
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'  Modlf  Icacions  to  improve  the  accuracy 
"f  the  skin  friction  and  the  shaoe 
parameter  of  the  velocity  profile  at 
,.'>w  '.-'cal  leynoids  numOer. 

(i)  Allowance  for  the  effect  on  turbu¬ 
lence  structure  of  flow  curvature  (a 
correction  which  Is  ivaiiabie  In  the 
lag-entrainment  method  but  not 
Implemented  In  vr,K'. 

^,1.1  Higher-order  effects 

jlfferenclng  the  respective  equations 
of  continuity  and  Reynolds-averaged, 
Navier-otokes  equations  for  two  flows  - 
tne  SqulvalenC  Inviscld  Flow  'HIP)  and  the 
Real  Viscous  Flow  (RV?")  -  East''*  was  able 
to  adapt  the  equations  for  Integrated  mass 
and  momentum  deficit  In  the  shear  layers 
to  Include  higher-order  effects.  The 
former  flow  coincides  with  the  latter 
wnere  the  flow  is  Inviscld  and  Is  the 
smooth  continuation  of  the  inviscld  flow 
wi‘:hin  the  shear  layers.  He  obtained  for 
the  wall  transpiration  velocity  In  EIP 

where  suffixes  I  and  w  refer  to  EIP 
and  to  wall  (aerofoil  contour  or  rear 
livldipg  streamline)  conditions,  fJ  la 
veLoclty  In  the  streamwise  direction  s 
and 

,  -5 

**  — 7] —  (2) 

oiw^^lw  ‘q  ^  ^ 

Is  ffene-aiiaed  displacement  thickness. 

The  Integration  Indicated  In  equation  (2) 
is  performed  from  the  wall  (z  *  0)  to  the 
outer  edge  of  the  shear  layer  Cz  *  5) 
along  a  normal  to  the  wall. 

Defining  a  generalised  momentum 
thickness 

1  ^ 

1  »  - 4-  /[pU(Ui„  -  IJ)  *• 

oiw'Jlw  0 

oiUlCUl  -  Idz 


Her9  Cf  Is  4n  equivalent  skin  mctlon 
coefficient,  M  Is  Mach  cumber  and 
H  •  l./e  .  The  terra  f  may  be  expressed 
as  the  sum  of  two  terms  due  respectively 
to  flow  curvature  and  to  Reynolds  normal 
stresses . 

Equations  ■:  1  ?  and  differ  from  the 
standard  first-order  expressions  In  vr,K  in 
that  the  Invlscld-flow  quantities  are 
defined  by  wail  conditions  In  EIE  rather 
than  by  those  at  the  outer  edi5e  of  the 
shear  layer  (which  In  vqK  are  Inferred 
from  the  wall  static  pressure  of  EVE  using 
the  boundary-layer  approx Ima t Ion 
3p/3z  •  0,  where  p  Is  static  pressure'. 
Lock  and  Firmin'^  argued  that  this  change 
Is  equivalent  to  an  Increase  In  adverse 
pressure-gradient  In  regions  of  rapid 
boundary-layer  growth  such  as  near  the 
trailing  edge.  An  additional  difference 
Is  that  the  right-hand  side  of 
equation  (3)  Is  non-zero.  In  general, 
whereas  It  Is  zero  In  the  first-order 
equation. 

The  second  term  on  the  right-hand  side 
of  equation  (3)  Is  generally  less  Import¬ 
ant  than  the  first  and  Is  thus  neglected. 
Consistent  with  this  approximation, 

Is  assumed  equal  to  the  surface  speed  in 
EIP,  q]^„  •  v^iw  +  '■'^iw  ^  study  of  the 

errors  In  boundary-lay“r  thicknesses 
Involved  In  using  Insteal  of 

suggests  that  they  may  be  Ignored  for  the 
flows  considered  In  this  paper. 

In  a  similar  way.  East  derived  the  nor¬ 
mal  momentum-integral  equation  and  from  It 
showed  that  wall  static  pressures  In  the 
two  flows  are  related  by  the  expression 

Plw  -  Pw  *  ‘*flw'^lw'9  *  «•)  . 

where  c»  Is  the  displacement-surf aoe 
curvature.  For  convenience,  this  equation 
la  not  used  directly  In  SVGK  but  Is  com¬ 
bined  with  the  usual  linearising  assump¬ 
tions  to  derive  a  flow  speed  corresponding 
to  p„ 


and  Ignoring  terms  which  are  significant 
only  when  surface  curvature  Is  large  and 
the  boundary  layer  grows  rapidly.  East 
showed  that  the  streamwlse  momentum 
Integral  equation  may  be  expressed  as 
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q„  -  U^^[l  +  x»(9  +  6»)1  .  (4) 

In  principle,  equation  (4)  allows  a 
solution  to  be  obtained  for  the  Inviscld 
flow  which  Is  compatible  with  the  Kutta 
condition  of  smooth  flow  at  the  trailing 
edge  expressed  as 

''’w^u  *  '’’wh  - 

suffixes  u  and  1  referring  to  the 
upper  and  lower  surfaces  approaching  the 
trailing  edge.  However,  the  conformal 
transformation  used  In  the  calculation  of 
the  Inviscld  flow  is  such  that  the  sol¬ 
ution  for  the  flow  Is  not  defined  at  the 
trailing  edge.  The  method  used  to  solve 


-.his  prohL*»m  within  the  fr^imework  of  the 
revi.'^e'l  treatment  of  the  i^hear  layers  Is 
leacrib*“d  later. 


'ieflnltlon  but  the  formula  used  ilffers 
from  that  of  t.he  lag-entralnment  method 
in  two  respects: 


Modlf Irac tons  to  the  entrainment  and  lag 
equations  In  the  lag-entrainment  method 
consequent  on  the  use  of  EIP  wail  quan¬ 
tities  In  the  definition  of  the  integral 
thicknesses  ar®  described  in  Ref  9- 

ohape  Parameter  Relationship 

In  the  lag  entrainment  method,  the  mass- 
flow  snao^  parameter 

0 

and  the  transformed  shape  parameter 

M  ^  MU,^  -  U)dz/p,^U^^8 
0 

are  related  by  the  eTpression 

‘I,  -  M5  ^  1.72/(H  -  1)  -  O.OKH  -  1)^  . 

This  formula  Is  Intended  to  represent  a 
wide  range  of  aerofoil  flows  but  Is  biased 
towards  those  with  severe,  adverse 
pressur”  gradients.  The  last  term  on  the 
right-hand  side  Is  Intended  to  ensure  that 
dPi^/dlT  Is  finite  for  all  positive  TT  , 
thus  avoiding  singular  behaviour  In  the 
■ntegratlon  of  the  shear-layer  equations 
by  the  direct  method. 

In  2VGK  an  alternative  relationship  Is 
used 

Hj  -  2  <■  l.!i/(H  -  1)  +  0.5(tT  -  1). 

1.3  <  fr  «  U  ; 

=  '1  t  (H  -  »),  U  <  H  <  12  .  (5) 


Equation  (5)  displays  a  minimum  In  H. 
at  a  value  of  H  (•  2.7)  close  to  that* 
for  Incipient  separation.  As  Implied 
above,  relationships  of  this  type  cannot 
be  used  for  the  calculation  of  separated 
flows  by  the  direct  method  normally 
associated  with  the  lag-entralnment  method 
but  pose  no  problem  when  the  Inverse 
method  Is  used.  LocK'*  showed  that  an 
expression  similar  to  equation  (5)  Is 
suitable  for  low-speed  flows,  and  It  would 
appear  that  the  same  Is  true  of  flows  at 
hlith  subsonic  speed  except  In  regions  of 
sudden  and  severe  pressure  gradient^''  (eg 
strong  shock  waves). 

A . 1 . 3  Skln-Frlctlon  Relationship 


fa)  Allowance  Is  made  for  the  effects  of 
low  Reynolds  number  on  the  velocity 
profile  of  the  turbulent  shear  layer. 

(b)  A  lower  limit  of  -0.0002  Is  Imposed 
on  the  skin  friction  coefficient. 

The  first  modification,  whlcn  Is 
described  In  detail  In  Ref  9.  Is  made  to 
allow  for  changes  In  the  character  of  tur¬ 
bulent  boundary-layers  at  low  Reynolds 
number  (ie  at  a  momentum-thickness 
Reynolds  number  R,  below  about  5000), 
relative  to  that  at  higher  Reynolds 
number,  originally  observed  by  Coles'’. 

According  to  Preston*’,  a  flat-plate 
boundary  layer  with  fully-developed  tur¬ 
bulence  is  not  possible  for  values  of  R, 
less  than  320  In  low  speed  flows.  In  view 
of  this,  and  In  the  absence  of  evidence  on 
the  effects  of  pressure  gradient  and 
compressibility,  this  value  has  been  taken 
to  be  the  lower  limit  of  R^  Just 
downstream  of  transition,  the  momentum 
thickness  there  being  adjusted  accordingly 
when  necessary. 

A.l.ii  Curvature  Effects  on  Turbulence 

Structure 


A  method  allowing  for  the  Influence  on 
turbulence  structure  of  streamwlse  flow 
curvature  Is  described  by  Sreen  et  ai > ’ . 
This  correction  with  a  modification  pro- 
posed  by  BradshawZ’  to  allow  for  a  lag  of 
10«  can  be  Included  In  3VCK.  The  mean 
value  of  flow  curvature  across  the  shear 
layer  used  In  this  correction  Is  taken  to 
be  that  of  the  displacement  surface.  It 
is  doubtful  if  this  approximation  Is  valid 
for  separated  flows  where  the  variation  of 
flow  curvature  across  the  layer  Is  likely 
to  be  large.  This  aspect  of  the  method  Is 
Considered  again  In  section  5  where  calcu¬ 
lation  and  measurement  are  compared. 

4.2  Vlsoous-Invlscld  Interaction 
Procedure 

BVGK  Is  seral-lnverse  In  character  since 
It  uses  an  Inverse  method  to  calculate  the 
turbulent  shear  layers  but  applies  the 
wall  transpiration  condition  In  the  calcu¬ 
lation  of  EIP  as  In  the  conventional 
direct  scheme.  In  both  VGK  and  BVGK  the 
method  of  Garabedlan  and  Korn  Is  used  to 
calculate  the  Invlscld  flow  with  a  modifi¬ 
cation  to  the  numerical  difference  scheme 
by  Lock*  to  Improve  the  representation  of 
shocks.  However,  In  BVGK  a  different 
method  Is  'ised  t;  satisfy  the  Kutta 
condition  and  this  Is  described  later. 


The  skln-frlctlon  coefficient  Cp  Is 
assumed  to  be  equivalent  to  the  standard 
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■1.2.1  Calculation  of  Shear  Layers  and 

Vlseoua-Invlacld  Matching 

’n  “aoh  surface,  the  Laminar  boundary 
layer  and  the  turbulent  boundary  layer  to 
one  i;rld  point  downstream  of  transition 
are  calculated  Ln  the  usual  direct  way. 
='urther  lownsr.resm,  the  Inverse  method  Is 
used  whereby  s*  Is  specified  and 
is  determined  along  with  other  boundary- 
layer  parameters,  the  Iteration  process 
being  started  with  an  assumed  distribution 
>f  dl.splacement  thlogness. 

'latching  between  the  Inverse  part  of  the 
solution  and  EIF  is  accomplished  with  an 
expression  due  to  Carter^  i 

x»<n  i-  U  ,  X»(h)  u  , 


where  superscripts  I  and  7  refer  to 
the  previous  Invlscld  30l".tlon  and  to  the 
solution  of  the  shear  layers,  n  Is  the 
number  of  viscous  Iterations  and  (1  Is  a 
relaxation  factor.  Wall  transpiration 
velocity  is  then  obtained  from 
epuaclon  (1)  with  oj.w  ^dd  being 

derived  from  the  inverse  solution  of  the 
shear  layers. 

11.2.2  Treatment  of  Trailing  Edge  and 

Wane 

The  conformal  mapping  used  in  the  calcu¬ 
lation  of  EIF  is  singular  at  the  trailing 
edge,  and  thus  the  solution  of  EIF  is  not 
defined  at  this  point,  ilenoa,  without 
further  consideration,  it  is  not  possible 
either  to  match  the  two  flow  solutions 
there  using  Carter's  formula  or  to  satisfy 
the  Kutta  condition  directly.  The  method 
used  to  solve  this  problem  Is  as  follows. 
First,  a  further  approximation  to  s*  at 
the  lralll.d(j  ’dgj  is  rp.j-H  pv  a  smooth 
interpolation  of  the  displacement  surface 
I  Fig  3).  For  this  purpose,  the  wake 
centre  line  near  the  trailing  edge  Is 
taken  to  be  a  straight  line,  Initially 
assumed  to  be  parallel  to  the  reference 
streamwise  axis  of  the  aerofoil.  Second, 
having  determined  displacement-surface 
curvature  at  the  Irailirg  edg-  on  each 
surface  by  Interpolation,  flow  speeds 
corresponding  to  RVP  pressures  at  the 
trailing  edge  are  found  using 
eduatlon  <>>).  At  this  stage,  the  Xutta 
condition  is  not  necessarily  satisfied; 
therefore  the  angle  of  the  wake  centre¬ 
line  relative  to  the  reference  axis,  *  , 
Is  adjusted  after  each  shear-layer  calcu¬ 
lation  until  values  of  5*  at  the 
trailing  edge  are  obtained  which  are  con¬ 
sistent  both  with  the  Xutta  condition  and 
the  requirement  of  a  smooth  displacement 
surf  ace  . 


•  $moatr«  fu/iaer 

•  HuMA  confiitiftn 


fig  3  Flow  conditions  to  be  satisfied  st 
trailing  edge 

The  condition  of  zero  pressure  differ¬ 
ence  across  the  wake  centre-line  is  also 
imposed  further  downstream.  As  shown 
by  Lock'^,  this  condition  defines  the  Jump 
in  velocity  across  the  wake  in  EIF. 


riq.4  Strfiewixe  rarnrion  of  jfteor  -  Ityer 
Ihiclincxs  near  friilinq  edqe  of  oluol- 
hot-  awrofoil 

The  effect  of  a  blunt  base  is  incorpor¬ 
ated  in  a  simple  way  as  follows.  As 
before,  the  displacement  surface  Is  taken 
to  be  continuous  at  the  trailing  edge,  and 
the  unique  shape-parameter  relationship, 
equation  (5),  Is  retained.  Pig  b 
Illustrates  the  streamwise  distributions 
of  the  shear-layer  parameters  upstream  and 
downstream  of  the  tra'Ilrg  idtr®  Fnth  the 
thickness  and  displacement  thickness 

of  the  shear  layer  above  the  rear 

dividing  streamline  Increase  dlsoontlnu- 
ously  at  the  trailing  edge  by  the  amount 
h  ,  the  height  of  the  trailing  edge  of  the 
upper  surface  above  the  Intersection  of 
the  dividing  streamline  with  the  base. 

This  Implies  that  Jumps  at  the 

trailing  edge  by  the  amount  h/e,j  ,  while 
(Hi)y  •  (1,,  -  Is  continuous  there. 

A  similar  argument  applies  to  the  shear 
layer  below  the  dividing  streamline 
{Pig  U).  Since  the  value  of  of 


7h‘=‘nr  Jowns r^a/n  if  ’:?5e 

“  ri  1 1  '  nj?  -'5  noc  neoessarlLy  ?on- 

'■'n.'v*:  '.'npi-e'i  '^y  -^nape- 

:Mr 'it; 'r  r“  I  t  ^  1  •: ns n  I p ,  an  art  J'la  trnenc  Is 
"  •)  ;{i  wn^ar®  necessary  as 

L  ;  ;js  t  n In  "Ig  The  position  :)f  the 

llvi'jirjg  atreamiln®  on  "he  hasej  «<hlcn 
1e".er'nlni“*5  "he  values  if  the  boundary- 
i.3v®r  ^hi ’'■<nesnes  either  5lde  the  >#ake, 

l3  inkno^yn;  however,  overall  forces  are 
fiund  to  ''e  Insensitive  to  ••’hangea  In  this 
position.  rhereCtre  In  the  oalculatlons 
'o  "hp  "'Lunt-oase  aerofoil  9AE  52i'J  tne 
oelght  )  !s  "OKen  to  oe  '}0%  tf  base 
t o  1  'knesn , 


■design)  Drrtlnates  were  used  In  "he  'al'. 
latlons;  the  errors  In  the  ordinates  an® 
small  with  a  'sagnltude  if  less  tnan 
■7.0002c,  where  :  la  aerofoi:.  ?nor'i.  ^h*- 
associated  errors  In  surfacp  pressur®  ar*’ 
calculated  tj  te  negligible  ■^T'^ept  :.osp 
to  sonic  condl"lon3. 

l)rag  Is  cal'.-ulated  jslng  "h®  expransl'" 
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■••(}  ;  4(Jigslwiffif  'a  04S5  '  Mow  1fH0« 
vs^  lawnjfream  of 
•flilinq  fdqf 

T-e  v<nK“  houndary  conditions  In  SI?  are 
ipnL!.'!  n  “.hs  cnordwlss  sxCsnslons  of  tbs 
• 'f  chn  bass  In  the  way  described  by 
’•’W  and  Will  lams  ^  ^  . 

3 _ lomparlson  Between  Caloulatlon 

and  Experiment 


?'Jll  letalls  of  the  calculation  pr-'- 
redur"  ar”  i5l'/en  In  ?ef  9.  Nominal 


Is  the  '/Iscous  drai?  coefficient,  suffix 
referring  to  conditions  far  lownstream. 

The  wave-dn^g  coefficient  'n,^  Is 

calculated  by  a  procedure^’  which  llent- 
Ifles  conditions  Just  upstream  of  the 
shock  In  the  flcwfleld  and  Infers  a  snocw 
total-pressure  loss  from  that  of  a 
Ranklne-Hugonlot  shock  of  the  same  Mach 
number  normal  to  and  just  upstream  of  the 
shook.  Wave  drag  then  follows  from  appli¬ 
cation  of  the  momentum  '.f.eorem  to  str“am- 
tubes  downstream  of  the  shocw  on  the 
assumption  of  adiabatic,  Isentroplc  flow. 

5.1  Convex  Sections 

Calculated  and  measured  pressure  Jlstrl- 
butlons  are  shown  In  ?lgs  b  and  '  for 
”  0.T3'',  C^  »  1.6  and  for  two  Reynolds 
numbers  R  •  sT)  »  lo'’,'or  I'.”  •  10°  for 
RAE  5230)  and  5  «  10'’.  These  flows  have  a 
supercrl t leal  -eglon  above  the  upper  sur¬ 
face  commencing  close  to  the  leading  edge 
and  terminating  In  a  weak  snook  at  about 
50t  chord.  In  some  cases,  re-expanslon  of 
the  flow  is  evident  upstream  of  the  shock. 
All  the  flows  are  sensitive  to  errors  In 


f'a6  Caltulatsd  and  measured  pressure  distributions  tor  Convex  sections.  ”.^.-0  135 


niz  "'y  '•/“r  ‘:ho3e  rf  VTK  13 

■' V  f:  ’:hy  .  F'^r  flows  without 

-M«  pr^a'1 1  1 'ns  ‘^C  irig 

"V  ^V'lK  i5‘?n*?ril  ly  wit;hln  ^bouC 

.’'<''0.  ■'f  -p«»  vsii'ii?.  ')r\  *:he 

"n«r  •’.an'i,  ic  -ilgper  :'r  lower  Incidences 
■1  n  ■' : r  mock-free  flow,  ohe  estl- 

?-es  M'  Py  ^V<1K  "are  not:  is  good, 

i '  0  '  r®pan  '  I -\*:  ilgr.'^r  Incidences 
>  '.n  ire  *^3  Dec  li  1 1  y  “7l'i<“nt  for  One 
.:w  "be  I irges f ,  regl  :'n  >r  leparaClon 

^A£  =  n  -»  10'^'.  As  ^o^^ed  oefore, 

/‘IK  “  InaPen  “.h»^  rear  loading  for 

icM  ■''-■'wn,  Trvi  nence  predl^ci  higner  suc- 
Irn.i  “.nin  of  rieasuremenC  opspream 

f  “n®  sn:  ck  in  i  given  lift  with  Phe  con- 
rno-^nce  “niP  "he  nerhod  overestimates 
av®  3rig.  "ome  Improvement  In  agreement 

=>^w*»en  'i.^'iLitlon  ind  measurement  Is 
")"ilnpd  f'-r  Ilf"  roef 1<:  lents  greater 
Mn  ino'it  ■'.*■'  ‘'or  TAE  ?230  it  R  =  6  ■«  10^ 
“he  ourvarjr®  correction  to  turbulence 
’  neglected  as  shown  In  ?lg  S. 


upper  Surface.  However,  3VGK  gives 
Improved  estimates  jf  pressure  llvStr 
butlons  compared  wltn  "hose  -it*  VdK, 
tlouiarly  at  "he  lower  rf  the  the  tw 
Reynolds  numbers  Fig  11'.  Figs  a  i 
show  that  3VGK  estimates  Irag  are 
significantly  closer  to  the  measured 
values  than  those  of  VTK.  Further, 
detailed,  comparisons  of  overall  for 
are  shown  In  Fig  11  for  "he  nominal 
number  and  Reynolds  numoer  >f  “Igs  ^ 
and  10.  Fig  11  tonflrms  the  ImDr'^ve 
accuracy  of  the  predictions  ^-f  Irag 
BVOK  compared  with  those  :f  7H?l .  -o 
sections,  3VGK  predictions  tf  drag  : 
riclenc  are  within  about  A .  MO  I  :f  " 
measured  values  over  the  range 
0.2  <  Of  <  0.6.  BVGK  Is  also  seen  t 


•  3VGK 

•  /GK 

s  =  S  .  -0*  .  . 


\r;'‘''ther  possible  source  of  error  arises 
'rcm  the  ss'^umptlon  of  the  wave-drag  pro- 
•edur®  "hat  the  local  flow  Is  normal  to 
-1,^  -n-'ck.  This  assumption  Is  llKeiv  to 
•‘*ru  t  Iv  an  overestimate  of  wave  drag  tn 
•'T''ws  wnere  the  shock  Is  oblique 
"h*>  aerofoil  surface. 
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Fressur*  1 1 s t r 1 bu t Ions  at  corresponding 
'rndMlons  "o  those  of  Figs  6  and  F  are 
In  Fig*'  9  and  10  for  the  relaxing 
-actions,  "he  flows  over  which  are  all 
rr-dl'to.}  cy  gvOK  to  be  attached  on  both 
5ur''ac«s.  RVGK  is  c^en  to  give  accurate 
®3tlmat‘»s  of  the  pressure  ils  tr  Ibtit  Ions 
aft  about  60t  cnord  but  discrepancies 
''etwe‘=*n  pr«dic’'*ons  by  QVGK  and  measure¬ 
ment  are  apparent  in  and  Just  downst.-eam 
tf  the  supercrltlcai-flow  region  on  the 
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“.he  Tiore  hohiirate  estimates  of  lift  at  a 
jlven  angle  of  Incidence  and  pltehlng- 
moment  at  a  gl'/en  lift. 

Comparison  Setween  Pigs  8  and  11  shows 
that  RAE  5236  has  larger  drag  at  a  given 
lift  chan  chat  of  the  convex  section 
RAE  5236.  This  occurs  In  spite  of  the 
flow  heing  calculated  to  be  attached  on 
the  former  section  and  separated  at  98* 
chord  on  the  upper  surface  of  Che  latter. 

5-3  Two-.i’art  3ectlon 

.“ressure  distributions  'for  M_  >  0.73, 

Cl.  '  0.6)  and  overall  forces  and  pitching 
moment  for  the  two-part  section  RAE  5234 
are  shown  In  Figs  12  and  13.  As  with  the 
other  sections,  BVOK  Is  seen  to  give 
accurate  estimates  of  drag  over  a  range  of 
lift  coefficients  which  Is  useful  for 
design  purposes  and  to  provide  signifi¬ 
cantly  Improved  predictions  of  pressure 
distributions,  lift,  drag  and  pitching 
moment  compared  with  those  of  VQK . 

Overall,  the  agreement  between  BVGK  pre¬ 
dictions  and  measurement  Is  good  for  all 
sections,  except  at  high  and  low  lift  when 
strong  shocks  appear.  Attention  Is  drawn, 
In  particular,  to  the  accuracy  of  the 
“Stlmates  of  differences  In  '  subcrltlcal ' 
drag  between  the  sections  of  all  three 
families  and  between  Reynolds  numbers  for 
a  given  section.  This  Illustrates  the 
possible  use  of  the  method  as  a  design 
tool  and  also  for  'extrapolating'  wind- 
tunnel  data  for  wing  sections  to  'full 
scale'  . 

6  Conclusions 

A  combined  theoretical  and  experimental 
Investigation  of  transonic  flows  over 
three  families  of  advanced  aerofoils  with 


varloift  rear-pressure  distributions  has 
been  described.  In  the  experiment, 
special  care  was  taken  to  ensure  that  the 
data  were  of  sufficient  quality  to  allow 
rigorous  validation  of  CPD  methods.  A 
large  number  of  different  flows  were 
studied,  ranging  from  those  that  were 
completely  attached  to  those  with  regions 
of  separation. 

The  data  have  been  used  to  validate  CFD 
methods  developed  at  RAE  including  a 
vlscous-lnvlscld  Interaction  method  known 
as  BVQK.  This  method  has  been  shown  to 
predict  accurately  pressure  distributions 
and  section  drag  for  a  wide  range  of  flows 
Including  those  with  rear  separation, 
suggesting  the  use  of  the  method  as  a  tool 
for  designing  sections  and  extrapolating 
wind-tunnel  data  to  'full  scale'. 

Following  extensive  validation,  BVOK  is 
currently  in  use  In  'JK  aerospace  Industry. 
Eventually,  the  method  will  be  replaced  by 
more  sophisticated  methods  but  before  they 
are  accepted  as  design  tools  they  will 
have  to  demonstrate  an  accuracy  at  least 
equal  to  that  shown  by  BVCK. 

References 

1  P..R.  Oarabedlan,  D.G.  Xorn.  "Analysis 
of  transonic  airfoils."  Comm.  Pure 
App  Math,  XXI'/,  pp  841-351,  (1971). 

2  M.R.  Collyer,  R.C.  lock.  "Prediction 
of  viscous  effects  In  steady  transonic 
flow  past  an  aerofoil."  utro  Qu,  30 
485,  (1979). 

3  R.H.  Doe,  A.  Pagano,  T.W.  Brown.  "The 
development  of  practical  Euler  .methods 
for  aerodynamic  design."  ICAS  86  5.5, 
(1936). 


13 


3.T.  'lien,  .  LI,  'I.  'Ihen. 

"''.liTi  iar.  Ion  of  transonic  airfoil 
flows  Oy  intaracClon  of  Suler  and 
oo'indary-iayer  equatlona.  '  AIAA  Paper 
>lo  37-  ,  {  1987  )  . 

L..J.  King.  "A  comparison  of  turbu¬ 
lence  oloaure  models  for  transonic 
flows  about  airfoils."  AIAA  Paper 
'Jo  87-0A13,  :  1987)  . 

L.-7.  Johnston.  "Some  preliminary 
re.su  Its  from  a  prediction  method  for 
the  viscous  flow  around  aerofoil 
sectlon.s."  ARA  Memorandum  38 1, 
i  10871  . 

T.L.  Holst.  "Viscous  transonic  air¬ 
foil  workshop  compendium  of  results." 
AIAA  Paper  No  87-l<i60  ,  (1987). 

L.F.  East,  P.D.  .Smith,  P.J.  Merryman. 
"Prediction  of  the  development  of 
separated  turbulent  boundary  layers  by 
the  lag-entralnment  method."  RAE 
Technical  Report  770JJ6,  (1977). 

P.R.  Ashlll,  R.P.  Wood,  D.J.  Weeks. 

"An  Improved,  semi-inverse  version  of 
the  viscous,  Sarabedlan  and  Korn 
metiiod."  RAE  Technical  Report  87003, 
n87)  . 

P.R.  Ashlll,  J.L.  Fulker,  D.J.  Weeks. 
"The  alr-lnjectlon  method  of  fixing 
boundary-layer  transition  and  Investi¬ 
gating  scale  effects."  The 
Aeronautical  .Journal  2ii  PP  3lJt-33tt, 

{  1987)  . 

P.R.  Ashlll,  D..J.  Weeks,  J.L.  Fulker. 
"Wind  tunnel  experiments  on  aerofoil 
models  for  the  assessment  of  compu¬ 
tational  flow  method."  AOARD-CPP- 
il37,  (  1988). 

P.R.  Ashlll,  D.J.  Weeks.  "A  method 
for  determining  wall-interference 
corrections  In  solid  wall  tunnels  from 
measurements  of  static  pressure  at  the 
wall."  AGARD-CP-335  Paper  1,  (1983). 

J.E.  Green,  D.J.  Weeks, 

J.W.F.  Brooman.  "Prediction  of  turbu¬ 
lent  boundary  layers  and  wakes  In  com¬ 
pressible  flow  by  a  lag  entrainment 
method."  ARC  RiM  3791,  (1973)- 

L.F.  East.  "A  representation  of 
second-order  boundary  layer  effects  In 
the  momentum  Integral  equation  and  In 
vlscous-lnvlscld  interactions."  RAE 
Technical  Report  81003,  (1981). 

,1.0.  Lock,  M.C.P.  Flrmln.  "Surveys  of 
techniques  for  estimating  viscous 
effects  In  external  aerodynamics.”  In 
"Numerical  methods  In  aeronautical 
fluid  dynamics"  (ed  P.L.  Roe), 


16  R.t:.  lock.  "Prediction  of  the  drag  of 
wings  at  subsonic  speeds  by 
vlscous/lnvlscld  interaction 
techniques."  AQARD-R-733  "Aircraft 
drag  prediction  and  reduction", 

{1985). 

17  D.J.  Weeks.  RAE  unpublished  work. 

18  D.E.  Coles.  "The  turbulent  boundary 
layer  In  a  compressible  fluid."  JSAF 
Project  Rand  Report  R-aOJ-PR,  '1963). 

19  J.H.  Preston.  "Minimum  Reynolds 
number  for  a  turbulent  boundary  layer 
and  the  selection  of  a  transition 
device."  Journal  of  Fluid  Mechanics 
3,  Ft  IV,  pp  373-38",  (1957)  . 

30  P.  Bradshaw.  "Effects  of  streamline 
curvature  on  turbulent  flow."  AGARD 

AG  169,  (1973)  . 

31  J.E.  Carter.  "A  new  boundary  layer 
Invlscld  Iteration  technique  for 
separated  flow."  AIAA  Paper  ’9-l{J50, 
(1979). 

33  R.C.  Look,  B.R.  Williams.  "Vlscous- 
Invlscld  Interactions  In  external 
aerodynamics."  Prog  Aerospace  3cl 
Vol  3<i,  pp  51-17)  ,  (  19871  . 

33  C.M.  Billing,  A.J.  Bocol.  "The 

MACHCONT  method  for  calculating  the 
wave  drag  of  a  3D  aerofoil."  ARA 
Memorandum  373,  (1986). 

Acknowledgment 

The  efforts  of  Mrs  N.  Rycrofl  and 

Mr  O.L.  Riddle  In  preparing  the  diagrams 

and  Miss  C.J.  Betts  In  assisting  with  "the 

calculations  are  gratefully  acknowledged. 


)M^,CoA»rQrOi)[l^aatot  CM(taU«r‘(Coiptt«<«.Aatliix>N|^aiid^l^^ 

7A7a00QB  _  Baytl  Atraipwai'  .  IK 


S*.  ^wtwii^  Atiney'tCode  6c  Spemotii%  A|«iK;  (6Mtti#tt  A«rthoifty)  NaiW  i^ 


7.  TW»  CAleulccioA  «Bd  pwurimint  b£  tranw^c  flews  bvsr  ssrofolls  with 
aavsl  rsto  saetienM 


7a.  fForTiaiiAatia«)  Titk  in  Fotatp  Isinune 


Th.  (For  Conference  Fapara)  Title,  Place  and  Date  of  Coofettnee 


f.Antbot  I.Sucnaiiw.Initiati  9a.  Authors 

Ashill*  t.  K. 


It.  Contract  Number 


12.  Period 


9b.  A«lliors3,4. 


13.  Prelect 


10.  Date  .  Pacet  Ref 

I  ’3  I23 


14.  Other  Reference  Nos. 


IS.  DWrteitianfUtanent 

(a) Cieiitroledby- 

(b)  afeciaHlniaatinna(if  aiy)  - 

if  RlcintMM  that  acoFjf  of  thhdoenaaaiU  AhD  be  lelcMed  oretieai  tefK  to  RAB  Leaflet  No3  to  SvMileneiit  6  of 

MOceiiiMcdd. 


16,  Deieripiort  (Ktywoida)  (DeKrlpton  marked  *  ate  eelectid  from  TEST) 

tMMosls.  Asrofolls 


IT.Abittaci 

A  irpiibiMiid  tbsvrscLesl  sad  sxpsiriMatAl  inowotisstioa  into  craaoeaie  flows 
o««*  MMoits  «£  adpsaeod  dMiga  is  daso^ad.  Aa  s^  haws  bsso  psrfoxaad 


seft«Nac  'S|*nads  sad  owar  a  srUa  saaga  «f  layaplds  utthtie  vg  to  20  aillioa  oa 
r  wicaAeils  aich  raat  o*’***«<^  9f  direar^  Joxa  sad 

stadtsd,  all  of  of  IdX 

if-giid-.iiiOa.fa  high  dagaaa' Of  >0*0.  fiiftisdrr  Soe  aaiofoti^' gniotna 

■IfOiMi.  'iiijid:o«ataU'::^f««a'«ta  OWgiiRld  .wjft  gtgdiBtloaa  h]r'a"’«aX«ttlat'ioB  ' 
igoiid.-oil;;  iM  oiaadjM“iadiogid  tdeodp^loa  ooooiltt-'i|^-:'iiwl^  al&oaaiaea  ia 

fa* '  - 


‘tmUMmM:. 


'lit^-giigadid 


